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The bilocal expansion of Borel transform provides an efficient way of Borel resummation with low order per- 
turbations in QCD. Its applications to the heavy quark pole mass, static potential, and lattice calculation are 
reviewed. 



1. Introduction 

The perturbative expansions in weak coupling 
constant in quantum field theories are in general 
divergent, with the coefficients growing in facto- 
rial. One way to make a rigorous sense out of the 
divergent series is to resum the series using Borel 
transform. When the large order behavior of the 
expansion is sign-alternating Borel resummation 
may yield the full amplitude. But in QCD the 
series are in general of same-sign and not Borel 
resummable, and must be augmented by nonper- 
turbative effects. However, Borel resummation 
can still be useful in these cases when a precise 
definition of the nonperturbative effects are re- 
quired, or the coupling is large enough that the 
usual perturbative amplitude from the low order 
perturbations is not sufficiently accurate for the 
purpose. 

The same-sign large order behavior of weak 
coupling expansion gives a singularity (renor- 
malon) on the contour of the Borel integral. 
When the coupling is small the Borel integral re- 
ceives most of its contribution from the immedi- 
ate neighborhood of the origin, where the Borel 
transform can be well-described by the low order 
perturbations, and the singular behavior of the 
Borel transform is not of much concern, but as 
the coupling increases the contribution from the 
region near the singularity becomes important, 
and one needs to have an accurate description of 
the Borel transform in the region that contains 
the origin as well as the nearest singularity to the 
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origin. 

It is obvious that it is a formidable task to re- 
build the singular behavior of the Borel trans- 
form using the ordinary perturbative expansions 
about the origin. In fact, the solvable models like 
the quantum mechanical instanton models sug- 
gest that one need many 10s orders of pertur- 
bative terms to obtain a good Borel resumma- 
tion. This approach is obviously not viable in 
QCD considering the high cost of loop calcula- 
tions. Thus an alternative approach is required 
to do any attempt for realistic Borel resumma- 
tion with the low order perturbations in QCD. 

The nature of the renormalon singularity can 
be determined by renormalization group tech- 
nique, and the (singular part of) Borel trans- 
form near the singularity can be systematically 
expanded perturbatively, up to the residue of the 
singularity that determines the overall normaliza- 
tion constant of the large order behavior (For a 
review see pQ). It was later realized that, once the 
nature of the singularity is known, the residue too 
can be calculated perturbatively, using only the 
usual perturbative expansions of the quantity in 
concern 0. 

We thus have a systematic expansion of the 
Borel transform at two points, the origin as well 
as the singularity. The bilocal expansion is to 
combine these expansions via an interpolation to 
obtain an accurate description of the Borel trans- 
form not only about the origin but also about the 
singularity 3 . 

The bilocal expansion has been successfully ap- 
plied to several QCD observables |3I4I5| . and fol- 
lowing is a review of those applications. 
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2. Bilocal expansion 

Let us assume that the Borel transform A(b) 
of an amplitude has the perturbative expansion 
about the origin to order M as 
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Table 1 

Residue for the Adler function with Nf quark flavor. 





Nf = 1 


Nf = 2 


N f = 3 


N f = 4 


LO 


0.904 


0.904 


0.946 


1.018 


NLO 


0.521 


0.546 


0.592 


0.674 


NNLO 


0.592 


0.549 


0.494 


0.411 



and expansion about the nearest renormalon sin- 
gularity to order N as 
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where bo is the location of the singularity. The 
"Analytic part" denotes terms that are analytic 
on the disk | fo — &o I < bo, and in general is not 
known. 

We can now combine the expansions JQ) and 
as a bilocal expansion using the following in- 
terpolation 0: 
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When we compare this expression to (J2J we see 
that the first sum simulates the "Analytic part" . 
The coefficients hi can be determined by demand- 
ing that the bilocal expansion, when expanded 
about the origin, reproduce the perturbative ex- 
pansion 0J. This yields the relation 

N 

h = a -j\f(l +2jcj) , 

3=0 



hi = a\ - —[1 + v + S^Cj(l + v - j)] 
oo , 
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The exact Borel transform is given by the limit 
of the bilocal expansion: 



A{b) = lim AuM b ) 

M,N^(X) 



(5) 



For the bilocal expansion to be usable it is 
necessary to compute the residue Af. Since the 
residue determines the normalization constant of 
the large order behavior, the Borel transform in 
bilocal expansion allows one to resum the di- 
vergent series to all orders once the residue is 
known. 

The residue can be perturbatively computed as 
follows. Noticing that 



M - R(bo) 
where 

R(b) = A(b)(l 



b/bo) u 



(G) 
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one can expand R(b) about the origin to evaluate 
JV. Notice that the expansion involves only the 
usual perturbative coefficients a*. Although R(b) 
is not analytic at bo, the expansion is nevertheless 
convergent because R(b) is bounded at the singu- 
larity. Fortunately, for several QCD observables 
this expansion gives a rather quick convergence. 
Adler function is such an example (see Table [TJ 

H- 

3. Pole mass 

The perturbative coefficients for the pole mass 
grow rapidly, and the convergence for the bottom 
quark, for example, is not good, as can be seen in 



^polc 



771^(1 + 0.093 + 0.045 + 0.032) (8) 



at a s = 0.22 , Nf = 4. 

The Borel transform for the heavy-quark pole 
mass has the nearest renormalon singularity at 
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6=1/2 and the bilocal expansion can be com- 
puted up to the order M = 2, N = 2. The pertur- 
bative calculation of the residue shows an excel- 
lent convergence 3 6 



M = 0.4244 + 0.1224 + 0.0101 = 0.5569 . 



(9) 



With this residue and N = 2, and varying M 
from to 2 we obtain the Borel resummed bottom 
quark pole mass ttibr as 

tobrMms = 1 + °- 1577 + °- 0041 ~ 0.0003 . (10) 

Notice the remarkable convergence. The leading 
term already contains the bulk of the contribu- 
tion, which shows that the proper handling of 
the singularity in the bilocal expansion acceler- 
ates the convergence in an enormous way. 

As the above calculation shows the Borel re- 
summation allows one to have a precise relation 
between the pole mass and MS mass. This means 
that the pole mass defined through Borel resum- 
mation is as physical as the MS mass and can 
be used as an alternative to the short distance 
masses. 

4. Static potential 

The heavy quark static potential obtained in 
perturbation theory disagrees badly with the lat- 
tice simulation as shown in Fig^ even in short 
distances where perturbation is expected to work. 

The nature of the nearest renormalon singular- 
ity is the same as the pole mass and the Borel re- 
summation can be performed in a similar fashion 
0. The result is given in Fig|2 The resummed 
potential agrees remarkably well with the lattice 
computation. This is a concrete example that 
shows the efficiency of the bilocal expansion. 

5. Lattice determination of heavy-quark 

mass 

The matching relation for the heavy-quark 
mass (b-quark) in lattice heavy quark effective 
theory reads |Zj 



pole 



A, A = £{a) — 8m(a) , 



(11) 



where A denotes the rcnormalizcd binding energy, 
a is the lattice spacing, Mb is the B-meson mass, 
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Figure 1. The static potential at leading order 
(dotted), next-leading order (dashed), and next-next- 
leading order (solid). The data points denote lattice 
potential. 



8m is the mass shift in the static limit, and £ 
is the binding energy that is to be computed in 
lattice simulation. 

The matching relation suggests that with high 
order perturbative calculations of the pole mass 
and the mass shift an accurate determination of 
the MS mass be possible with precision calcula- 
tion of £. There is, however, a problem with that 
the pole mass and mass shift do not converge well 
perturbatively, due to the renormalon. To over- 
come this difficulty the matching relation is ex- 
pressed in terms of the MS mass only, bypassing 
the pole mass, as 



MS 



A(o)(l 



where A (a) = M B - £{a). 

Although this expression does not suffer from 
the renormalon divergence there could be a con- 
vergence problem coming from the mixing of the 
two independent scales, the lattice spacing and 
the heavy quark mass, through the RG scale \x. 
When the two scales are close, as in the bottom 
quark case, this problem may not be severe but 
when the two scales are far-separated this may 
become a serious problem. 

The scale mixing problem can be solved with 
Borel resummation. The matching relation lllj l 
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Figure 2. Borel resummed potential 



with Borel resummation reads 



0.8 



M R - A 



Bli 



(13) 



where A BR = £{a) — <5m BR (a). Here the pole 
mass and mass shift are resummed independently, 
using the strong coupling at a scale optimal for 
each quantities. The scales are cleanly isolated, 
with 1/a only in Sm BR and m\, only in m BR . 

Table |21 shows the MS mass for the bottom 
quark from the above two approaches in quenched 
limit. In this case the results from the two ap- 
proaches agree remarkably well. However, when 
the lattice spacing and the heavy-quark mass are 
far-different the scale mixing problem appears. 
For an (imaginary) heavy-light meson with mass 
Tfiff — 500 GeV we obtain, using (|13|) . the BR 
mass: 

roq R = 499.542 , 499.542 , 499.523 (GeV) (14) 
at 1/a = 2.12, 2.91, 3.85 (GeV), respectively, and 



Table_2 

The MS mass m™ s determined in the Borel summa- 
tion method and PCM. The units are in GeV. 



a- 1 


2.12 


2.91 


3.85 


BR method 


4.312 


4.312 


4.297 


PT method 


4.319 


4.320 


4.311 



using the relation 

to qV to cF= 1 + 0.03638 + 0.00008-0.00004 (15) 
we obtain the corresponding MS mass 



MS 



= 481.984 , 481.984 , 481.965 (GeV) . (16) 



What is remarkable with this result is that the 
uncertainties in these numbers are only ±20 MeV, 
entirely coming from A BR and the conversion to 
Wq S from the BR mass, both of which cause less 
than 20 MeV uncertainty. 

On the other hand, with the conventional 
method using (|12|) . the TOq S obtained shows 
strong dependence on the RG scale fi which can 
be anywhere in a wide range between 1/a and 
Mq, and the uncertaintities in the heavy quark 
mass from this approach are several hundred 
MeVs, which are an order of magnitude larger 
than those in the Borel resummation method. For 
details we refer the readers to 

6. Conclusion 

With the bilocal expansion Borel resummation 
for certain QCD observables can be done to a 
remarkable accuracy with the low order pertur- 
bations. The salient feature of the Borel resum- 
mation is the automatic isolation of hierarchical 
scales, which can be useful in quarkonium sys- 
tems. 
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